APS/123-QED 



First-principles calculation of the Gilbert damping parameter via the linear response 
formalism with application to magnetic transition-metals and alloys 
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A method for the calculations of the Gilbert damping parameter a is presented, which based on 
the linear response formalism, has been implemented within the fully relativistic Korringa-Kohn- 
Rostoker band structure method in combination with the coherent potential approximation alloy 
theory. To account for thermal displacements of atoms as a scattering mechanism, an alloy-analogy 
model is introduced. This allows the determination of a for various types of materials, such as 
elemental magnetic systems and ordered magnetic compounds at finite temperature, as well as for 
disordered magnetic alloys at T = K and above. The effects of spin-orbit coupling, chemical and 
temperature induced structural disorder are analyzed. Calculations have been performed for the 
3d transition-metals bcc Fe, hep Co, and fee Ni, their binary alloys bec Fei-iCoi, fee Nii-^Fea,, 
fee Nii-ajCoa; and bcc Fei- x V x , and for 5d impurities in transition-metal alloys. All results are in 
satisfying agreement with experiment. 

PACS numbers: 72.25.Rb 71.20.Bc 71.70.Ej 75.78.-n 



I. INTRODUCTION 

During the last decades dynamical magnetic properties 
have attracted a lot of interest due to their importance in 
the development of new devices for spintronics, in par- 
ticular, concerning their miniaturization and fast time 
scale applications. A distinctive property of such devices 
is the magnetization relaxation rate characterizing the 
time scale on which a system being deviated from the 
equilibrium returns to it, or how fast the device can be 
switched from one state to another. In the case of dy- 
namics of a uniform magnetization /vecM this property 
is associated with the Gilbert damping parameter G(M) 
used first in the phenomcnological Landau-Lifshitz (LL) 1 
and Landau-Lifshitz-Gilbert (LLG) theory 2 describing 
the magnetization dynamics processes by means of the 
equation: 



ldM 1 , 

— = -M x H off + M x 

7 dr 



G(M) dM 
7 2 M| ~dT 



(1) 



where M s is the saturation magnetization, 7 the gy- 
romagnctic ratio and H e // = — c?M.F[M(r)] being the 
effective magnetic field. Sometimes it is more conve- 
nient to use a dimensionless Gilbert damping parame- 
ter a given by a = G/{^M S ) (see, e.g. 3-5 ). Safonov 
has generalized the Landau-Lifshitz equation by intro- 
ducing a tensorial form for the Gilbert damping parame- 
ter with the diagonal terms characterising magnetization 
dissipation 6 . Being introduced as a phenomcnological pa- 
rameter, the Gilbert damping is normally deduced from 
experiment. In particular, it can be evaluated from the 
resonant line width in ferromagnetic-resonance (FMR) 
experiments. The difficulty of these measurements con- 
sists in the problem that there exist several different 
sources for the broadening of the line width, which have 



been discussed extensively in the literature 7-13 . The line 
width that is observed in ferromagnetic resonance spec- 
tra is usually caused by intrinsic and extrinsic relax- 
ation effects. The extrinsic contributions are a conse- 
quence of spatially fluctuating magnetic properties due 
to sample imperfections. Short range fluctuations lead 
to two magnon scattering while long range fluctuations 
lead to an inhomogencous line broadening due a super- 
position of local resonances 14 . In order to separate the 
intrinsic Gilbert damping from the extrinsic effects it is 
necessary to measure the frequency and angular depen- 
dence of the ferromagnetic resonance line width, e. g. 
two magnon scattering can be avoided when the mag- 
netization is aligned along the film normal 11 (perpen- 
dicular configuration). Usually one finds a linear fre- 
quency dependence with a zero frequency offset and one 
can write AH (10) = + AH(0). When such measure- 
ments are performed over a wide frequency range the 
slope of AH as a function of frequency can be used 
to extract the intrinsic Gilbert damping constant. In 
metallic ferromagncts Gilbert damping is mostly caused 
by electron magnon scattering. In addition Gilbert-like 
damping can be caused by eddy currents. The magni- 
tude of the eddy current damping is proportional to d 2 , 
where d is the sample thickness 10 . In sufficiently thin 
magnetic films (d < 10 nm) the eddy current damping 
can be neglected 10 . However, for very thin films relax- 
ation mechanisms that occur at the interfaces can also 
increase and even dominate the damping. Such effects 
are spin pumping 15,16 and the modified electronic struc- 
ture at the interfaces. In the present work spin pumping 
and the modified interface electronic structure are not 
considered and we assume that bulk-like Gilbert damp- 
ing dominates. 

Much understanding of dynamical magnetic properties 
could in principle be obtained from the simulation of 
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these processes utilizing time-dependent first-principles 
electronic structure calculations, that in turn would pave 
the way to developing and optimizing new materials for 
spintronic devices. In spite of the progress in the de- 
velopment of time-dependent density functional theory 
(TD-DFT) during the last decades 17 that allows to study 
various dynamical processes in atoms and molecules from 
first principles, applications to solids are rare. This is due 
to a lack of universally applicable approximations to the 
exchange-correlation kernel of TD-DFT for solids. Thus, 
at the moment, a tractable approach consists in the use of 
the classical LLG equations, and employing parameters 
calculated within a microscopic approach. Note however 
that this approach can fail dealing with ultrafast mag- 
netization dynamics, which is discussed, for instance, in 
Refs. [18 and 19], but is not considered in the present 
work. 

Most of the investigations on the magnetization dissi- 
pation have been carried out within model studies. Here 
one has to mention, in particular, the so-called s-d or 
p-d exchange model 20-23 based on a separate considera- 
tion of the localized 'magnetic' d-electrons and dclocal- 
ized s- and p-electrons mediating the exchange interac- 
tions between localized magnetic moments and responsi- 
ble for the magnetization dissipation in the system. As 
was pointed out by Skadsem et al. 24 , the dissipation pro- 
cess in this case can be treated as an energy pumping 
out of the d-electron subsystem into the s-electron bath 
followed by its dissipation via spin-flip scattering pro- 
cesses. This model gave a rather transparent qualita- 
tive picture for the magnetization relaxation in diluted 
alloys, e.g. magnetic semiconductors such as GaMnAs. 
However, it fails to give quantitative agreement with ex- 
periment in the case of itinerant metallic systems (e.g. 
3d- metal alloys), where the d-states are rather delocaliscd 
and strongly hybridized with the sp-electrons. As a con- 
sequence the treatment of all valence electrons on the 
same footing is needed, which leads to the requirement 
of first-principles calculations of the Gilbert damping go- 
ing beyond a model-based evaluation. 

Various such calculations of the Gilbert damping pa- 
rameter are already present in the literature. They usu- 
ally assume a certain dissipation mechanism, like Kam- 
bersky's breathing Fermi surface (BFS) 25,26 , or more gen- 
eral torque-correlation models (TCM) 3,27 . These models 
include explicitly the spin-orbit coupling (SOC), high- 
lighting its key role in the magnetization dissipation pro- 
cesses. However, the latter methods used for electronic 
structure calculations cannot take explicitly into account 
disorder in the system that in turn is responsible for the 
aforementioned spin-flip scattering process. Therefore, 
this has to be simulated by using external parameters 
characterizing the finite lifetime of the electronic states. 
This weak point was recently addressed by Brataas et al. A 
who described the Gilbert damping by means of scatter- 
ing theory. This development supplied the formal basis 
for the first parameter-free investigations on disordered 
alloys for which the dominant scattering mechanism is 



potential scattering caused by chemical disorder 5 . 

Theoretical investigations of the magnetization dissi- 
pation by means of first-principles calculations of the 
Gilbert damping parameter already brought much un- 
derstanding of the physical mechanisms responsible for 
this effect. First of all, key roles are played by two ef- 
fects: the SOC of the atomic species contained in the 
system and scattering on various imperfections, either 
impurities or structural defects, phonons, etc. Account- 
ing for the crucial role of scattering processes respon- 
sible for the energy dissipation, different types of scat- 
tering phenomena have to be considered. One can dis- 
tinguish between the ordered-compound or pure-element 
systems for which electron-phonon scattering is a very 
important mechanism for relaxation, and disordered al- 
loys with dominating scattering processes resulting from 
randomly distributed atoms of different types. In the 
first case, the Gilbert damping behavior is rather differ- 
ent at low and high temperatures. At high temperature 
atomic displacements create random potentials leading 
to SOC-induced spin-flip scattering. At low temperature, 
where the magnetization dissipation is well described via 
the BFS (Breathing Fermi-surface) mechanism 25 ' 26 , the 
spin-conserving electron-phonon scattering is required to 
bring the electronic subsystem to the equilibrium at ev- 
ery step of the magnetization rotation, i.e. to reoccupy 
the modified electronic states. 

In this contribution we describe a formalism for the 
calculation of the Gilbert damping equivalent to that 
of Brataas et oi. 4 , however, based on the linear rc- 
sponce theory 28 as implemented in fully rclativistic mul- 
tiple scattering based Korringa-Kohn-Rostokcr (KKR) 
formalism. It will be demonstrated that this allows to 
treat elegantly and efficiently the temperature depen- 
dence of a in pure crystals as well as disordered alloys. 



II. THEORETICAL APPROACH 

To have direct access to real materials and to obtain 
a deeper understanding of the origin of the properties 
observed experimentally, the phenomenological Gilbert 
damping parameter has to be treated on a microscopic 
level. This implies to deal with the electrons responsible 
for the energy dissipation in the magnetic dynamical pro- 
cesses. Thus, one equates the corresponding expressions 
for the dissipation rate obtained in the phenomenologi- 
cal and microscopic approaches E mag — E^is- Although 
a temporal variation of the magnetization is a required 
condition for the energy dissipation to occur, the Gilbert 
damping parameter is defined in the limit ui — > (see 
e.g., Ref. [24]) and therefore can be calculated within the 
adiabatic approximation. 

In the phenomenological LLG theory the time depen- 
dent magnetization M(i) is described by Eq. (1). Ac- 
cordingly, the time derivative of the magnetic energy is 
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given by: 



mag 



= H cff ~ = -l(rii) T [G(M)iii] 
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(2) 



where m = M/M s denotes the normalized magnetiza- 
tion. To represent the Gilbert damping parameter in 
terms of a microscopic theory, following Brataas et al. 4 , 
the energy dissipation is associated with the electronic 
subsystem. The dissipation rate upon the motion of the 

magnetization E c \ is = is determined by the under- 

lying Hamiltonian H(t). Assuming a small deviation of 
the magnetic moment from the equilibrium the normal- 
ized magnetization rii(r) can be expanded around the 
equilibrium magnetization rho 



m(r) = m + u(r) 



(3) 



resulting in the expression for the linearized time depen- 
dent Hamiltonian for the system brought out of equilib- 
rium: 

9 -H(mo) . (4) 



H = ^ (m ) + ^tv — 

UMi, 




Ei)S(E F - E. 



Due to the small deviation from the equilibrium, Ea\ s can 
be obtained within the linear response formalism, leading 
to the expression 4 : 

Edis = — Trfi ^ ^ u^iiv 

ij )iv 

xS(E F 

where Ep is the Fermi energy and the sums run over 
all eigenstates of the system. As Eq. (5) characterizes 
the rate of the energy dissipation upon transition of the 
system from the tilted state to the equilibrium, it can 
be identified with the corresponding phcnomenological 
quantity in Eq. (2), E mas = Edis- This leads to an ex- 
plicit expression for the Gilbert damping tensor G or 
equivalent ly for the damping parameter a = G/(^M S ) 
(Ref. [4]): 



fry 



dH 



du„ 




dH 



xS(E F — Ei)S(E F — Ej) , (6) 

where the summation is running over all states at the 
Fermi surface Ep. 

In full analogy to the problem of electric 
conductivity 29 , the sum over eigenstates \ipi) may be 
expressed in terms of the retarded single-particle Green's 
function lmG+{E F ) = -tt£. ^{^{Ep - E t ). This 
leads for the parameter a to a Kubo-Greenwood-likc 
equation: 




Trace 



dH 



Im G+(E F ) — Im G + (E F ) 



(7) 



with (...) c indicating a configurational average in case of 
a disordered system. 

The most reliable way to account for spin-orbit cou- 
pling as the source of Gilbert damping is to evaluate 
Eq. (7) using a fully relativistic Hamiltonian within the 
framework of local spin density formalism (LSDA) 30 : 



H = col ■ p + f3mc 2 + V(v) + (3cr ■ mJB(r) 



(8) 



Here ct^ and (5 are the standard Dirac matrices, a de- 
notes the vector of relativistic Pauli matrices, and p is 
the relativistic momentum operator 31 . The functions 
V(r) and B = cr • m_B(r) are the spin-averaged and 
spin-dependent parts, respectively, of the LSDA poten- 
tial. The spin density m s (r) as well as the effective ex- 
change field B(r) are assummed to be collinear within 
the unit cell and aligned along the z-direction in the 
equilibrium (i. e. m Sl o(r) = m s (r)mo = m s (r)e 2 and 
B (r) = £?(r)rh = B(r)e z ). Tilting of the magnetiza- 
tion direction by the angle 9 according to Eq. (3), i.e. 
m s (r) = m s (r)m = m s (r) (sin 9 cos <f>, sin 9 sin <j>, cos 9) 
and B(r) = _B(r)m leads to a perturbation term in the 
Hamiltonian 

AV(r) = /3tr ■ (m - m )B(r) = /3cr ■ uB(r) , (9) 
with (see Eq. (4)) 



■^-H(mo) = /3<7 M B(r) 
du a 



(10) 



The Green's function G + in Eq. (7) can be obtained in 
a very efficient way by using the spin-polarized relativis- 
tic version of multiple scattering theory 30 that allows us 
to treat magnetic solids: 



G+(r, r', E) = ^ Zftr, E) (E) Z™* (r', E) 

AA' 

S nm E t Z A(r, E) (r', E)Q(r' n - r n ) 

A 

+J2(r,E)ZZ?(r',E)e(r n -r' n )] (11) 

Here r, r' refer to site n and m, respectively, where 
Z£(r,E) = Z A (r n ,E) = Z A (r - R n ,E) is a function 
centered at site R„. The four-component wave functions 
Z^(r,E) (J^(r,E)) are regular (irregular) solutions to 
the single-site Dirac equation labeled by the combined 
quantum numbers A (A = (k, /x)), with k and /x being 
the spin-orbit and magnetic quantum numbers 31 . The 
superscript x indicates the left hand side solution of the 
Dirac equation, rj^™ (E) is the so-called scattering path 
operator that transfers an electronic wave coming in at 
site m into a wave going out from site n with all possible 
intermediate scattering events accounted for. 

Using matrix notation with respect to A, this leads to 
the following expression for the damping parameter: 

= -J— ^2 Trace (T° M r°™ f r"°) c (12) 

with the g- factor 2(1 + /x or b//x sp m) in terms of the spin 
and orbital moments, [i av i n and /J, rb, respectively, the 
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total magnetic moment \it 



f^spin ~l~ l^orbi 7"aA 



5?( T AA' — t A'a) an< ^ with tne ener gy argument E F omit- 
ted. The matrix elements in Eq. (12) are identical to 
those occurring in the context of exchange coupling 32 : 



T»? A = A2J x (r) 



d 



fl-(mo) 



(13) 



The expression in Eq. (12) for the Gilbert-damping 
parameter a is essentially equivalent to the one obtained 
within the torque correlation method (see e.g. Refs. [33- 
35]). However, in contrast to the conventional TCM the 
electronic structure is here represented using the retarded 
electronic Green function giving the present approach 
much more flexibility. In particular, it does not rely on 
a phenomenological relaxation time parameter. 

The expression Eq. (12) can be applied straightfor- 
wardly to disordered alloys. This can be done by de- 
scribing in a first step the underlying electronic struc- 
ture (for T = OK) on the basis of the coherent po- 
tential approximation (CPA) alloy theory. In the next 
step the configurational average in Eq. (12) is taken fol- 
lowing the scheme worked out by Butler 29 when dealing 
with the electrical conductivity at T = K or residual 
resistivity, respectively, of disordered alloys. This im- 
plies in particular that so-called vertex corrections of the 
type (T^mG+T,MiG+) c - (T fl lmG+) c (T,lmG+) c that 
account for scattering-in processes in the language of the 
Boltzmann transport formalism are properly accounted 
for. 

One has to note that the factor in Eq. (12) is sep- 
arated from the configurational average (•••)„, although 
both values, g and fJ-tot, have to represent the average per 
unit cell doing the calculations for compounds and al- 
loys. This approximation is rather reasonable in the case 
of compounds or alloys where the properties of the ele- 
ments of the system are similar (e.g. 3d-element alloys), 
but can be questionable in the case of systems containing 
elements exhibiting significant differences (3d-5d-, 3d-4f- 
compounds, etc), or in the case of non- uniform systems 
as discussed by Nibarger et al 36 . 

Thermal vibrations as a source of electron scattering 
can in principle be accounted for by a generalization of 
Eqs. (7) - (13) to finite temperatures and by including 
the electron-phonon self-energy E e ;_ p /j when calculating 
the Green's function G + . Here we restrict our considera- 
tion to elastic scattering processes by using a quasi-static 
representation of the thermal displacements of the atoms 
from their equilibrium positions. The atom displaced 
from the equilibrium position in the lattice results in a 
corresponding variation At™ = t n — ig 1 of the single-site 
scattering matrix in the global frame of reference 37,38 . A 
single-site scattering matrix t n (the underline denotes a 
matrix in an angular momentum representation A) for 
the atom n displaced by the value s™ from the equilib- 
rium position in the lattice can be obtained using the 



transformation matrices 37,39 



Ul LI {s v ,E)=^Y, il " +l ~ 1 ' 



X C WL nj l ,,( S yE)Y L n(^) . (14) 

Here m e is the electron mass, ji a spherical Bessel func- 
tion, Cll'L" stands for the Gaunt coefficients, and a 
non-relativistic angular momentum representation with 
L = (I, mi) has been used. Performing a Clcbsch-Gordon 
transformation for the transformation matrix U™ L , to 
the relativistic A representation, the t matrix t™ for the 
shifted atom can be obtained from the non-shifted one 
tg from the expression 



(15) 



Treating for a discrete set of displacements s™ each 
displacement as an alloy component, we introduce an 
alloy-analogy model to average over the set s" that 
is chosen to reproduce the thermal root mean square 
average displacement \J (u 2 )t for a given temperature 
T. This in turn may be set according to (u 2 )t = 
\^ez\^tP + i] with ^d/T) the Debye func- 
tion, h the Planck constant, k the Boltzmann constant 
and o the Debye temperature 40 . Ignoring the zero tem- 
perature term 1/4 and assuming a frozen potential for 
the atoms, the situation can be dealt with in full analogy 
to the treatment of disordered alloys on the basis of the 
CPA (see above). 

For small displacements the transformation Eq. (14) 
can be expanded with respect to s ra (see Ref. [39]) re- 
sulting in a linear dependence on s™ for non-vanishing 
contributions with Al = \l — l'\ = ±1. This leads, in 
particular, in the presence of atomic displacements for 
transition- metals (TM), for which an angular momen- 
tum cut-off of Z max = 2 in the KKR multiple scattering 
expansion is in general sufficient for an undistorted lat- 
tice, to an angular momentum expansion up to at least 
= 3. However, this is correct only under the assump- 
tion of very small displacements allowing linearisation of 
the transformation U with respect to the displacement 
amplitude s. Thus, since the temperature increase leads 
to a monotonous increase of s, the cut-off l max should also 
be increased. 



III. MODEL CALCULATIONS 

In the following we present results of calculations for 
which single parameters have artificially been manipu- 
lated in the first-principles calculations in order to sys- 
tematically reveal their role for the Gilbert-damping. 
This approach is used to disentangle competing influ- 
ences on the Gilbert-damping parameter. 
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A. Vertex corrections 

The impact of vertex corrections is shown in Fig. 1 
for two different cases: Fig. 1(a) represents the Gilbert 
damping parameter for an Fc±- x V x disordered alloy as a 
function of concentration, while Fig. 1(b) gives the cor- 
responding value for pure Fe in the presence of temper- 
ature induced disorder and plotted as a function of tem- 
perature. Both figures show results calculated with and 
without vertex corrections allowing for comparison. First 
of all, a significant effect of the vertex corrections is no- 
ticeable in both cases, although the variation depends on 
increasing concentration of V in the binary Fci-a; Va; alloy 
and the temperature in the case of pure Fe, respectively. 
Some differences in their behavior can be explained by 
the differences of the systems under consideration. Deal- 
ing with temperature effects via the alloy analogy model, 
the system is considered as an effective alloy consisting 
of a fixed number of components characterizing different 
types of displacements. Thus, in this case the tempera- 
ture effect is associated with the increase of disorder in 
the system caused only by the increase of the displace- 
ment amplitude, or, in other words - with the strength 
of scattering potential experienced by the electrons rep- 
resented by t n (T) — Vq. In the case of a random alloy 
the A\- X B X variation of the scattering potential, as well 
as the difference V\ — tfg, upon changing the concentra- 
tions is less pronounced for small amounts of impurities 
B and the concentration dependence is determined by the 
amount of scatterers of different types. However, when 
the concentration of impurities increases, the potentials 
of the components are also modified (this is reflected, 
e.g., in the shift of electronic states with respect to the 
Fermi level, that will be discussed below) and this can 
lead to a change of the concentration dependence of the 
vertex corrections. An important issue which one has to 
stress that neglect of the vertex corrections may lead to 
the unphysical result, a < 0, as is shown in Fig. 1(a). 
In terms of the Boltzmann transport formalism, this is 
because of the neglect of the scattcring-in term 41 lead- 
ing obviously to an incomplete description of the energy 
transfer processes. 



B. Influence of spin-orbit coupling 

As was already discussed above, the spin-orbit coupling 
for the electrons of the atoms composing the system is 
the main driving force for the magnetization relaxation, 
resulting in the energy transfer from the magnetic subsys- 
tem to the crystal lattice. Thus, the Gilbert damping pa- 
rameter should approach zero upon decreasing the SOC 
in the system. Fig. 2 shows the results for Py+15%Os, 
where \fa is plotted as a function of the scaling param- 
eter of the spin-orbit coupling 42 applied to all atoms in 
the alloy. As one can see, y/a has a nearly linear depen- 
dence on SOC implying that a varies in second order in 
the strength of the spin-orbit coupling 43 . 
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FIG. 1. The Gilbert damping parameter for (a) bec Fei- x ~V x 
(T — K) as a function of V concentration and (b) for bcc-Fe 
as a function of temperature. Full (open) symbols give results 
with (without) the vertex corrections. 
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FIG. 2. The Gilbert damping parameter for Py+15%Os as 
a function of the scaling parameter of spin-orbit coupling ap- 
plied to all atoms contained in the alloy. Red dashed line in 
plot - linear fit. The values and 1 for the SOC scaling pa- 
rameter correspond to the scalar-relativistic Schrodinger-like 
and fully relativistic Dirac equations, respectively. 
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IV. RESULTS AND DISCUSSIONS 
A. 3d transition-metals 

We have mentioned above the crucial role of scatter- 
ing processes for the energy dissipation in magnetiza- 
tion dynamic processes. In pure metals, in the absence 
of any impurity, the electron-phonon scattering mecha- 
nism is of great importance, although it plays a different 
role in the low- and high-tempcraturc regimes. This was 
demonstrated by Ebert et al. 2S using the alloy analogy 
approach, as well as by Liu 44 et al. using the 'frozen ther- 
mal lattice disorder' approach. In fact both approaches 
are based on the quasi-static treatment of thermal dis- 
placements. However, while the average is taken by the 
CPA within the alloy analogy model the latter requires 
a sequence of super-cell calculations for this purpose. 

As a first example bec Fe is considered here. The cal- 
culations have been performed accounting for the tem- 
perature induced atomic displacements from their equi- 
librium positions, according to the alloy analogy scheme 
described in section II. This leads, even for pure systems, 
to a scattering process and in this way to a finite value 
for a (see Fig. 3(a)). One can see that the experimen- 
tal results available in the literature are rather different, 
depending on the conditions of the experiment. In par- 
ticular, the experimental results Expt. 2 (Rcf. [45]) and 
Expt. 3 (Rcf. [46]) correspond to bulk while the measure- 
ments Expt. 1 (Rcf. [47]) have been done for an ultrathin 
film with 2.3 nm thickness. The Gilbert damping con- 
stant obtained within the present calculations for bec Fe 
(circles, a = 5.44 a.u.) is compared in Fig. 3(a) with 
the experiment exhibiting rather good agreement at the 
temperature above 100 K despite a certain underestima- 
tion. One can also sec a rather pronounced increase of 
the Gilbert damping observed in the experiment above 
400 K (Fig. 3(a), Expt. 2 and Expt. 3), while the theo- 
retical value shows only little temperature dependent be- 
havior. Nevertheless, the increase of the Gilbert damp- 
ing with temperature becomes more pronounced when 
the temperature induced lattice expansion is taken into 
account, that can be seen from the results obtained for 
a = 5.45 a.u. (squares). Thick lines are used to stress 
the temperature regions for which corresponding lattice 
parameters are more appropriate. At low temperatures, 
below 100 K, the calculated Gilbert damping parameter 
goes up when the temperature decreases, that was ob- 
served only in the recent experiment 47 . This behavior is 
commonly denoted as a transition from low-temperature 
conductivity-like to high-tempcraturc resistivity-like be- 
havior reflecting the dominance of intra- and intcr-band 
transitions, respectively 3 . The latter are related to the 
increase of the smearing of electron energy bands caused 
by the increase of scattering events with temperature. 
Note that even a small amount of impurities reduces 
strongly the conductivity- like behavior 28,45 , leading to 
the more pronounced effect of impurity-scattering pro- 
cesses due to the increase of scattering events caused by 



chemical disorder. Large discrepancies between the lat- 
ter experimental data 47 and theoretical results of the a 
calculations for bec Fe are related to the very small thick- 
ness of the film investigated experimentally, that leads to 
an increase of spin-transfer channels for magnetization 
dissipation as was discussed above. 

Results for the temperature dependent Gilbert- 
damping parameter a for hep Co are presented in Fig. 
3(b) which shows, despite certain underestimation, a rea- 
sonable agreement with the experimental results 45 . The 
general trends at low and high temperatures are similar 
to those seen in Fe. 

The results for pure Ni are given in Fig. 3(c) that show 
in full accordance with experiment a rapid decrease of a 
with increasing temperature until a regime with a weak 
variation of a with T is reached. 

Note that in the discussions above we have treated a 
as a scalar instead of a tensorial quantity ignoring a pos- 
sible anisotropy of the damping processes. This approx- 
imation is reasonable for the systems considered above 
with the magnetization directions along a three- or four- 
fold symmetry axis (see, e.g., the discussions in Rcf. [48 
and 49]). For a more detailed discussion of this issue 
the anisotropy of the Gilbert damping tensor a(M) has 
been investigated for bec Fe. To demonstrate the depen- 
dence of a on the magnetization direction M, the cal- 
culations have been performed for M = £|M| with the 
z axis taken along the (001), (111) and (011) crystallo- 
graphic directions. Fig. 4 presents the temperature de- 
pendence of the diagonal elements a xx and a yy . As to be 
expected for symmetry reasons, a xx differs from a yy only 
in the case of zj|(011). One can see that the anisotropic 
behavior of the Gilbert damping is pronounced at low 
temperatures. With an increase of the temperature the 
anisotropy nearly disappears, because of the smearing of 
the energy bands caused by thermal vibrations 49 . A sim- 
ilar behavior is caused by impurities with a random dis- 
tribution, as was observed for example for the Feo.95Sio.05 
alloy system. The calculations of the diagonal elements 
a xx and a yy for two different magnetization directions 
along (001) and (011) axes give a xx = a yy = 0.00123 in 
the first case and a xx = 0.00123 and a yy = 0.00127 in 
the second, i.e. the damping is nearly isotropic. 

The damping parameter a increases very rapidly with 
decreasing temperature in the low temperature regime 
(T < 100 K) for all pure ferromagnetic 3d metals, Fe, 
Co, and Ni (see Fig. 3), leading to a significant discrep- 
ancy between theoretical and experimental results in this 
regime. The observed discrepancy between theory and 
experiment can be related to the exact limit u> = taken 
in the expression for the Gilbert damping parameter. Ko- 
renmann and Prange have analyzed the magnon damp- 
ing in the limit of small wave vector of magnons q — > 0, 
assuming indirect transitions in the electron subsystem 
and taking into account the finite lifetime r of the Bloch 
states due to electron-phonon scattering. They discuss 
the limiting cases of low and high temperatures showing 
the analogy of the present problem with the problem of 
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FIG. 3. Temperature variation of the Gilbert damping pa- 
rameter of pure systems. Comparison of theoretical results 
with experiment: (a) bcc-Fe: circles and squares show the re- 
sults for ideal bcc Fe for two lattice parameters, a = 5.42 a.u. 
and a — 5.45 a.u.; stars show theoretical results for bcc Fe 
(a = 5.42 a.u.) with 0.1% of vacancies (Expt. 1 - Ref. [47], 
Expt. 2 - Ref. [45], Expt. 3 - Ref. [46]); (b) hep-Co: circles 
show theoretical results for ideal hep Co, stars - for Co with 
0.03% of vacancies, and 'pluses' - for Co with 0.1% of vacan- 
cies (Expt. Ref. [45]); and (c) fcc-Ni (Expt. Ref. [45]). 



extreme cases for the conductivity leading to the normal 
and anomalous skin effect. On the basis of their result, 
the authors point out that the expression for the Gilbert 
damping obtained by Kambersky 25 , with a ~ r is cor- 
rect in the limit of small lifetime (i.e. qvFT <C 1, in their 
model consideration, where q is a magnon wave vector 
and «f is a Fermi velocity of the electron). In the low- 
temperature limit the lifetime r increases with decreas- 
ing T and one has to use the expression corresponding 
to the 'anomalous' skin effect for the conductivity, i.e. 
a ~ tan -1 (qvpr) / qvp , leading to a saturation of a upon 
the increase of t. 

Another possible reason for the low-temperature be- 
havior of the Gilbert damping observed experimentally 
can be structural defects present in the material. To 
simulate this effect, calculations have been performed for 
fee Ni and bcc Fe with 0.1% of vacancies and for hep 
Co with 0.1% and 0.03% of vacancies. Fig. 3(a)-(c) 
shows the corresponding temperature dependence of the 
Gilbert damping parameter approaching a finite value for 
T — » 0. The remaining difference in the T-dependent be- 
havior can be attributed to the non-linear dependence of 
the scattering cross section at low temperatures as is dis- 
cussed in the literature for transport properties of metals 
and is not accounted for within the present approxima- 
tion. 



B. 3d Transition-metal alloys 

As is mentioned above, the use of the linear response 
formalism within multiple scattering theory for the elec- 
tronic structure calculations allows us to perform the nec- 
essary configurational averaging in a very efficient way 
avoiding supercell calculations and to study with mod- 
erate effort the influence of varying alloy composition 
on a. The corresponding approach has been applied to 
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FIG. 4. Temperature variation of the a„ and a yy com- 
ponents of the Gilbert damping tensor of bcc Fe with the 
magnetization direction taken along different crystallographic 
directions: M = z|M| ]|(001) (circles), M||(011> (squares), 
M||(lll) (diamonds). 
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FIG. 5. (a) Theoretical results for the Gilbert damping pa- 
rameter of bcc Fer-iCoa; as a function of Co concentration: 
CPA results for the bcc structure (full circles) describing 
the random alloy system, results for the partially ordered 
system (opened square) for x — 0.5 (i.e. for Fei-ajCo^ 
alloy with CsCl structure and alloy components randomly 
distributed in two sublattices in the following proportions: 
(Feo.9Coo.i)(Feo.iCoo.<)), the NL-CPA results for random al- 
loy with bcc structure (opened circles) and the NL-CPA re- 
sults for the the system with short-range order within the 
first-neighbor shell (opened diamonds). The dashed line rep- 
resents the DOS at the Fermi energy, Ef, as a function of Co 
concentration, (b) spin resolved DOS for bcc Fei-^Coa for 
x = 0.01 (dashed line) and x = 0.5 (solid line). 



the ferromagnetic 3d-transition-mctal alloy systems bcc 
Fei-zCo^, fee Nii.^Fe^, fee Ni^Co^ and bcc lei , V, . 

Fig. 5(a) shows as an example results for the Gilbert 
clamping parameter a{x) calculated for bcc ¥e\- x Co x 
for T = K at different conditions. Full circles rep- 
resent the results of the single-cite CPA calculations 
characterizing the random Fc-Co alloy. These results 
are compared to those obtained employing the non-local 
CPA 52,53 (NL-CPA) assuming no short-range order in 
the system (opened circles). Dealing in both cases (CPA 
and NL-CPA), with completely disordered system, the 
NL-CPA maps the alloy problem on that of an impurity 
cluster embedded in a translational invariant effective 
medium determined selfconsistently, thereby accounting 
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FIG. 6. The Gilbert damping parameter for Fei-zCoz (a) 
Nii-^Co^ (b) and Nii-^Fez (c) as a function of Co and Fe 
concentration, respectively: present results within CPA (full 
circles), experimental data by Oogane 50 (full diamonds), (d) 
Results for bcc Fer-xVa as a function of V concentration: 
T = 0K (full circles) and T = 300 K (open circles). Squares: 
experimental data 51 . Open circles: theoretical results by 
Starikov et al. 5 . 
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for nonlocal correlations up to the range of the cluster 
size. The present calculations have been performed for 
the smallest NL-CPA clusters containing two sites for bcc 
based system, accounting for the short-range order in the 
first-neighbor shell. As one can see, this results in a small 
decrease of the a value in the region of concentrations 
around x = 0.5 (opened diamonds), that is in agreement 
with the results obtained for partially ordered system 
(opened square) for x = 0.5. The latter have been calcu- 
lated for the Fei_ x Co2; alloy having CsCl structure and 
alloy components randomly distributed in two sublattices 
in the following proportions: (Feo.gCbo.iXFeo.iCoo.g). 

Because the moments and spin-orbit coupling strength 
do not differ very much for Fe and Co, the variation of 
a(x) should be determined in the concentrated regime 
essentially by the electronic structure at the Fermi en- 
ergy Ep. As Fig. 5(a) shows, there is indeed a close 
correlation with the density of states n{Ep) that may be 
seen as a measure for the number of available relaxation 
channels. The change of a{x) due to the increase of the 
Co concentration is primarily determined by an appar- 
ent shift of the Fermi energy also varying with concen- 
tration (Fig. 5(b)). The alloy systems considered have 
the common feature that the concentration dependence 
of a is governed by the concentration dependent density 
of states n{Ep), 

A comparison of theoretical a values with the experi- 
ment for bcc Fei-^Co^ is shown in Fig. 6(a), demonstrat- 
ing satisfying agreement. In the case of Ni^^Fe^ and 
Ni 1 _ 2 ,Co 2 , alloys shown in Fig. 6, (b) and (c), the Gilbert 
damping decreases monotonously with the increase of the 
Fe and Co concentration, in line with experimental data. 
At all concentrations the experimental results are under- 
estimated by theory approximately by a factor of 2. The 
calculated damping parameter a(x) is found in very good 
agreement with the results based on the scattering theory 
approach 5 demonstrating numerically the equivalence of 
the two approaches. An indispensable requirement to 
achieve this agreement is to include the vertex correc- 
tions mentioned above. As suggested by Eq. (12) the 
variation of a(x) with concentration x may also reflect 
to some extent the variation of the average magnetic mo- 
ment of the alloy, fitot- 

The peculiarity of the Fei- x V x alloy when compared to 
those discussed above is that V is a non-magnetic metal 
and has only an induced spin magnetic moment. De- 
spite that, the concentration dependence of the Gilbert 
damping parameter at T = K for small amounts of 
V (see Fig. 6(d)) displays the same trend as the pre- 
viously discussed alloys shown in Fig. 6(a)-(c). Taking 
into account a finite temperature of T = 300 K changes 
a value significantly at small V concentrations leading 
to an improved agreement with experiment for pure Fe, 
while it still compares poorly with the experimental data 
at xy = 0.27. One should stress once more that the con- 
centration dependent behavior of the Gilbert damping 
parameter of the alloys discussed above is different for 
an increased amount of impurities (more than 10%), as a 



result of a different variation of the DOS n(Ep) caused by 
a concentration dependent modification of the electronic 
states and shift of the Fermi level. 



C. 5d impurities in 3d transition metals 

As discussed in our recent work 28 investigating the 
temperature dependent Gilbert damping parameter for 
pure Ni and for Ni with Cu impurities, a is primarily 
determined by the thermal displacement in the regime of 
small impurity concentrations. This behavior can also be 
seen in Fig. 7, where the results for Fe with 5d-impurities 
are shown. Solid lines represent results for T = K for 
an impurity concentration of 1% (full squares) and 5% 
(full circles). As one can see, at smaller concentrations 
the maximum of the Gilbert damping parameter occurs 
for Pt. With increasing impurity content the a parame- 
ter decreases in such a way that at the concentration of 
5% a maximum is observed for Os. 

The reason for this behavior lies in the rather weak 
scattering efficiency of Pt atoms due to a small DOS 
n(Ep) of the Pt states when compared for example for 
Os impurities (see Fig. 9). This results in a slow decrease 
of a at small Pt concentration when the BFS mechanism 
is mostly responsible for the energy dissipation. A con- 
sequence of this feature can be seen in the temperature 
dependence of a (T = 300 K, opened squares): a most 
pronounced temperature induced decrease of the a value 
is observed for Pt and Au. When the concentration of 
5d- impurities is increased up to 5%, the maximum in a 
occurs for the element with the most efficient scatter- 
ing potential resulting in spin-flip scattering processes 
responsible for dissipation. The temperature effect at 
this concentration is very small. 

Considering in more detail the temperature dependent 
behavior of the Gilbert damping parameter for Fe with 
Os and Pt impurities, shown in Fig. 8, one can also ob- 
serve the consequence of the features mentioned above. 
At 1% of Pt impurities a decreases much steeper upon 
increasing the temperature, as compared to the case of 
Os impurities. Therefore, in the first case the role of the 
scattering processes due to atomic displacements is much 
more pronounced than in the second case with rather 
strong scattering on the Os impurities. When the con- 
centration increases to 5% the dependence of a on the 
temperature in both cases becomes less pronounced. 

The previous results can be compared to the results for 
the 5d- impurities in the permalloy FC80N120 (Py), which 
has been investigated also experimentally 54 . This system 
shows some difference in the concentration dependence 
when compared with pure Fe, because Py is a disordered 
alloy with a finite value of the a parameter. Therefore, 
a substitution of 5d impurities leads to a nearly linear 
increase of the Gilbert damping with impurity content, 
just as seen in experiment 54 . 

The total damping for 10% of 5d-impuritics shown in 
Fig. 10(a) varies roughly parabolically over the 5d TM se- 
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FIG. 7. Gilbert damping parameter for bcc Fe with 1% 
(squares) and 5% (circles) of 5d impurities calculated for 
T = OK (full symbols) and for T = 300K (opened sysmbols). 



ries. This variation of a with the type of impurity corre- 
lates well with the density of states n 5d (EF) (Fig. 10(b)). 
Again the trend of the experimental data is well repro- 
duced by the calculated values that are however some- 
what too low. 



V. SUMMARY 

In summary, a formulation for the Gilbert damping pa- 
rameter a in terms of linear response theory was derived 
that led to a Kubo-Greenwood-like equation. The scheme 
was implemented using the fully rclativistic KKR band 
structure method in combination with the CPA alloy the- 
ory. This allows to account for various types of scat- 
tering mechanisms in a parameter-free way, that might 
be either due to chemical disorder or to temperature- 
induced structural disorder (i.e. electron-phonon scat- 
tering effect). The latter has been described by using 
the so-called alloy-analogy model with the thermal dis- 
placement of atoms dealt with in a quasi-static manner. 
Corresponding applications to pure metals (Fc, Co, Ni) 
as well as to disordered transition-metal alloys led to very 
good agreement with results based on the scattering the- 
ory approach of Brataas et al. and well reproduces the 
experimental results. The crucial role of vertex correc- 
tions for the Gilbert damping is demonstrated both in 
the case of chemical as well as structural disorder and 
the accuracy of finite-temperature results is analyzed via 
test calculations. 

Furthermore, the flexibility and numerical efficiency 
of the present scheme was demonstrated by a study 
on metallic systems on a scries of binary 3d-alloys 
(Fci-zCoz, Nii^Fcz, Nii-zCoz and Fei_ x V x ), 3d- 5d 
TM systems, the permalloy-5<i TM systems. The agree- 
ment between the present theoretical and experimental 
results is quite satisfying, although one has to stress 
a systematic underestimation of the Gilbert damping 
by the numerical results. This disagreement could be 
caused either by the idealized system considered theoret- 
ically (e.g., the boundary effects arc not accounted for 
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FIG. 10. (a) Gilbert damping parameter a for Py/5d TM 
systems with 10 % 5d TM content in comparison with 
experiment 54 ; (b) spin magnetic moment m^ in and density 
of states u(Ef) at the Fermi energy of the 5d component in 
Py/5d TM systems with 10 % 5d TM content. 
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